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ABSTRACT 


A  flutter  theory  based  on  modified  Yates  hydrodynamic 
loading  was  used  to  predict  the  flutter  characteristics  of 
five  hydrofoil  flutter  models.  Theoretical  flutter  speeds 
were  unconservative  when  flutter  was  predicted  to  occur  in 
a  predominantly  bending  mode,  and  conservative  when  flutter 
was  predicted  to  occur  in  a  predominantly  torsional  mode. 

The  theory  may  be  useful  for  hydrofoils  which  are  susceptible 
to  flutter  in  a  torsional  mode. 

ADMINISTRATIVE  information 

This  work  was  authorized  and  funded  under  the  Hydrofoil  Development 
Program  of  the  Naval  Ship  Systems  Command,  Subproject  S4606,  Task  1703. 

INTRODUCTION 

The  phenomemon  of  flutter  involves  mutual  interactions  among  hydro- 
dynamic,  elastic,  and  inertial  forces;  therefore,  a  flutter  potential 
exists  for  any  system  that  is  subjected  to  these  forces.  Hydrofoil  flutter 
is  conceptually  identical  to  aeroelastic  flutter  of  aircraft  lifting-sur¬ 
faces.  However,  hydrofoils  differ  significantly  from  airfoils  in  the 
magnitude  of  certain  parameters  relating  flow  characteristics  to  structural 

properties--mass  ratio  and  reduced  frequency.  The  magnitude  of  mass  ratio 

1  2 

is  a  determining  factor  in  the  accuracy  of  flutter  theory.  * 

For  typical  airplane  wings,  mass  ratio  is  of  the  order  of  50;  for 
hydrofoils  and  ship  control  surfaces,  mass  ratio  is  approximately  0.5.  In 
the  high  mass-ratio  region  typical  of  airfoils,  flutter  theories  employed 
in  the  aeronautical  field  in  general  have  proven  to  be  adequate.  When 
these  theories  were  extended  to  the  low  mass -ratio  region  typical  of  hydro¬ 
foils,  however,  flutter-speed  predictions  were  almost  invariably  uncon¬ 
servative.  Various  studies  have  led  to  the  belief  that  the  classical 

2  3  4 

formulations  of  unsteady  hydrodynamic  loads  are  not  adequate. *■’  ’ 

A  promising  approach  was  initiated  by  Yates ^  when  he  applied  a 
striplike  structural  and  hydrodynamic  load  theory  to  relatively  low 
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mass-ratio  foils.  Yates  used  a  load  theory  which  had  been  modified  to 
account  for  three-dimensional  flow  effects  by  inserting  calculated  span- 
wise  distributions  of  lift  slope  and  aerodynamic  center.  This  approach  had 
been  proven  successful  in  predicting  aerodynamic  flutter.6  When  applied  to 
a  hydrofoil  flutter  model  of  mass  ratio  0.99,  the  Southwest  Research  Insti¬ 
tute  (SwRI)  30-inch  model,4  the  Yates  theory  predicted  a  flutter  speed  which 
was  conservative  by  20  percent.  In  view  of  this  conservative  prediction  at 
a  low  mass  ratio,  as  well  as  other  favorable  results,  it  was  felt  that  the 
theory  should  be  investigated  further. 

In  his  discussion  of  results,  Yates*  suggested  that  the  remaining 
difference  between  theory  and  experiment  might  be  due  to  the  incomplete 
nature  of  his  loading  modification.  The  striplike  loading  expressions  had 
incorporated  three-dimensional  effects  in  the  circulatory  loading  terms; 
however,  the  noncirculatory  loading  terms  had  been  left  in  two-dimensional 
form.  This  formulation  had  proved  adequate  in  the  aerodynamic  mass-ratio 
range,  where  the  noncirculatory  terms  contributed  a  small  fraction  of  the 
total  loading.  However  the  noncirculatory  terms,  which  include  virtual 
mass,  become  relatively  more  important  as  mass  ratio  decreases.  Therefore 
a  three-dimensional  modification  to  the  noncirculatory  portion  of  the  load¬ 
ing  might  improve  the  accuracy  of  hydrofoil  flutter  predictions. 

Following  the  Yates  suggestions,  the  authors  further  modified  the 
unsteady  loading  expressions  used  by  Yates  to  incorporate  three-dimensional 
effects  in  the  noncirculatory  terms.  A  spanwise  modification  factor  was 
developed  from  two-dimensional  virtual  mass  expressions  and  was  applied  to 
the  noncirculatory  terms  to  form  a  fully  modified  unsteady  loading  theory. 
The  loading  was  then  combined  with  a  highly  versatile,  lumped-parameter 
structural  representation  to  form  a  theory  capable  of  predicting  dynamic 
response  and  flutter  characteristics  of  a  wide  range  of  winglike  structural 
configurations. 

To  evaluate  the  fully  modified  strip  theory,  1 .ve  different  flutter 
models  were  analysed  for  which  experimental  results  were  available.  The 
mass  ratios  of  the  models  ranged  from  0.09  to  0.99.  All  models  were 
cantilevered  hydrofoils;  four  of  the  five  were  surface-piercing,  while  the 
fifth  was  fully  submerged. 


The  analytical  treatment  of  each  model  consisted  of  calculating  the 

steady  loading  for  the  model,  using  three-dimensional  lifting  surface 
7 

theory;  calculating  the  noncirculatory  load  distribution,  using  the  expres¬ 
sion  developed  by  the  authors;  and  solving  the  general  flutter  equations. 

To  avoid  convergence  difficulties  and  numerical  instability,  which  had 

3  8  9 

presented  problems  in  other  investigations,  *  *  the  flutter  equations  were 
solved  directly  rather  than  by  the  usual  normal  mode  approach. 

This  report  describes  the  complete  theoretical  flutter  calculations, 
including  the  noncirculatory- load  modification  and  the  method  of  structural 
analysis.  The  five  flutter  models  are  described,  and  experimental  and 
theoretical  results  are  compared. 

FLUTTER  ANALYSIS 


FLUTTER  THEORY 

Hydrofoils  susceptible  to  flutter  are  generally  slender,  vinglike 
structures  of  high  aspect  ratio.  It  is  possible  to  represent  such  struc¬ 
tures  by  using  a  straight  elastic  axis  with  structural  properties  lumped 
at  discrete  points  along  the  axis.  Hydrodynamic  forces  can  be  similarly 
lumped  at  stations  along  the  axis. 

In  the  present  work,  each  hydrofoil  model  treated  was  divided  into 
10  sections.  Finite-difference  equations  were  written  for  the  structural 
and  fluid  interactions  among  the  foil  sections  according  to  the  procedure 
shown  in  Figure  1  and  described  in  the  next  paragraph.  Calculations  not 
presented  in  this  report  showed  that  the  differences  in  solutions  between 
a  foil  divided  into  10  sections  and  one  divided  into  20  sections  were 
negligible  for  the  first  three  modes,  implying  that  a  10-section  representa¬ 
tion  gave  a  converged  solution. 

The  governing  equations  for  flutter  analysis  are  derived  by  applying 
D'Alembert's  Principle  at  each  lumped-parameter  station  with  proper  con¬ 
straint  relations  to  account  for  boundary  conditions.  Figure  2  shows  each 
lumped-parameter  station.  The  flutter  equations  may  be  written  in  matrix 
form  in  the  following  manner 


[M] 


\h.\  (h.j 

[h . 

4  !C1  Kj 

+  (1  +  jg)  [K] 

!  i 

l9J 

(1) 


3 


where  is  the  independent  coordinate  to  represent  the  bending  degree  of 
freedom  at  the  ith  lumped-mass  station,  and 

is  the  independent  coordinate  to  represent  the  torsional  degree  of 
freedom  at  the  ith  lumped-mass  station. 

The  left  side  of  Equation  (1)  represents  the  structural  character¬ 
istics  of  the  hydrofoil.  Matrix  elements  for  the  mass,  damping,  and  stiff¬ 
ness  matrices  associated  with  beamlike  structures  such  as  the  present 
hydrofoils  may  be  derived  as  described  in  Reference  10. 

The  right  side  of  Equation  (1)  represents  the  hydrodynamic  force 
applied  to  the  hydrofoil.  The  force  may  be  expressed  in  the  same  form  as 
the  structural  representation 


*  [C*] 


(2) 


Substituting  Equation  (2)  into  Equation  (1),  one  may  arrive  at  the 
following  governing  matrix  equation  to  represent  the  hydroelastic  system 


[Ml 


♦  [C] 


=  {0} 


where  [M]  =  [M]  +  [M*] 

[C]  =  [C]  +  [C*] 

[K]  =  (1  +  jg)  [K]  *  [K*] 


(3) 


By  assuming 


h. 

l 


Hiest,  and  0.  =  0.est 


and  substituting  them  into  Equation  (3),  one  may  arrive  at  the  following 


expressions 


♦  s[C]  +  [K] 


=  {0} 


(4) 


This  represents  a  typical  complex  eigenvalue  and  eigenvector  problem. 
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Solving  Equation  (4),  one  may  obtain  the  complex  eigenvalues  of  s 
and  their  corresponding  eigenvectors  at  different  speeds.  Each  eigenvalue 
corresponding  to  a  mode  of  oscillation,  may  be  written  in  the  form 

s  =  Re(s)  +  jlm(s) 


or 


at  speed  V 


(5) 


where  Pe(s)-is  ths'real  part  of  s  or  the  negative  of  the  exponential  decay 
factor  *  -<$,  and 

Im(s)  is  the  imaginary  part  of  s  or  the  damped  natural  frequency. 

The  variation  of  Re(s)  of  all  the  modes  as  a  function  of  speed  may 
be  obtained  from  the  previously  described  calculation.  Flutter  is  con¬ 
sidered  to  occur  at  the  lowest  speed  for  which  Re(s)  becomes  zero.  The 
corresponding  flutter  frequency  and  flutter  mode  are  given  by  Im(s)  and 
its  eigenvector,  respectively. 

The  dynamic  characteristics  of  the  structure  in  vacuum  and  in  water 
at  zero  speed  may  be  obtained  from  Equation  (3)  by  omitting  certain  terms. 
First,  by  neglecting  all  the  unsteady  hydrodynamic  forces  and  all  damping 
elements.  Equation  (3)  becomes 


Solving  Equation  (6) ,  one  may  obtain  the  natural  frequencies  of  the  struc¬ 
ture  and  its  associated  mode  shapes  in  vacuum.  Second,  by  including  the 
virtual  mass  matrix  [M*]  in  Equation  (6),  one  obtains 


The  solution  of  Equation  (7)  gives  the  natural  frequencies  and  associated 
mode  shapes  of  the  structure  in  water  at  zero  speed. 

A  digital  computer  program**  has  been  used  to  solve  Equations  (4), 
(6),  and  (7)  on  the  CDC  6700  computer  at  the  Center. 
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FLUTTER  MODELS 


The  hydrofoils  selected  for  this  study  were  the  SwRI  30-inch  flutter 

4  8  9 

model,  Grumman  Models  1,2,  and  3,  and  Grumman  Model  A.  All  the  models 

had  very  low  mass  ratio,  ranging  from  0.09  to  0.99.  All  except  the  SwRI 

model  had  been  tested  as  partly  submerged  foils  perpendicular  to  the  water 

surface.  The  SwRI  model  had  been  tested  fully  submerged  with  a  reflection 

plate  at  its  root.  The  detailed  characteristics  of  each  model  are  summarized 

in  the  Appendix. 

HYDRODYNAMIC  LOADING 

As  shown  in  Equation  (2) ,  expressions  for  hydrodynamic  lift  and 
moment  are  required  as  functions  of  the  foil  displacements  h  and  0  and  the 
corresponding  velocities  and  accelerations.  In  the  original  Yates  treat¬ 
ment,  the  circulatory  part  of  the  classical  Theodorsen  two-dimensional 
12 

load  theory  was  modified  to  include  calculated  or  measured  spanwise  lift- 
slope  and  aerodynamic -center  distributions.  In  accordance  with  the  Yates 
suggestion,5  the  authors  made  an  additional  modification  to  the  loading 
used  by  Yates.  The  additional  modification  consisted  of  a  reduction  in  the 
magnitude  of  the  noncirculatory  loading  along  the  span  of  loaded  surfaces 
in  order  to  account  for  three-dimensional  flow  effects.  A  reduction  factor 
resembling  the  distribution  of  virtual  mass  on  a  flat  plate  was  chosen.  In 
the  following  sections,  the  location  of  the  modification  in  the  Yates  load¬ 
ing  expressions  is  given,  and  the  method  of  calculating  the  modification 
is  described.  In  addition,  the  method  is  described  for  calculating  the 
lift-slope  and  aerodynamic-center  distributions  required  for  the  circulatory 
modification,  and  the  calculated  distributions  are  given. 

MODIFIED  YATES  LOADING 

The  unsteady  hydrodynamic  expressions  which  Yates  applied  to  each 
strip  along  the  span  of  the  foil  were  similar  in  form  to  the  expressions 
presented  in  Reference  13  for  the  lift  and  pitching  moment  on  a  section  of 
a  wing  oscillating  harmonically  in  two-dimensional  incompressible  flow  at  a 
nonzero  sweep  angle.  Yates  reformulated5  the  expressions  to  permit  the  use 
of  arbitrary  values  of  steady  lift-curve  slope  and  section  aerodynamic  center, 
instead  of  the  values  2tt  and  quarter  chord,  respectively, 
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the  results  of  two-dimensional  thin  airfoil  theory.  The  required  spanwise 
distributions  of  section  lift-curve  slope  and  section  aerodynamic  center 
for  a  particular  planform  may  be  obtained  by  any  reliabile  theoretical  or 
experimental  method.  For  the  present  analysis  they  were  obtained  from  the 
kernel -function  type  of  lifting-surface  theory  presented  in  Reference  7. 

Modification  of  noncirculatory  loading  was  made  with  a  multiplicative 
factor  p^,  inserted  in  the  Yates  unsteady  loading  expressions  as  follows 


-P.  =  p.trpb.2  I’h.  -  V  0.  +  V  a.  tan  A  +  b.a.  (0.  ♦  V  T.  tan  A  j 
l  ri  l  li  ni  ni  ea  ii\i  ni  ea/ 


-C*  ,PVnb.  C(k)  w. 

a,i 


-M.  =  p.irpb.  vl/8  +  a.2j(0.  +  V  i.  tan  A  )  +  p.npb.2V  (h.  +  V  a.  tan  A  ) 
l  *i  \  i/\i  ni  ea/  *i  l  n\  l  n  l  ea/ 

+  d  Trpb.3a./h.  +  V  a. tan  A  J  -  p.irpb.2  V  2  (e.-a.b.T.  tan  A  ^ 
pi  M  i  i  \  l  n  i  ea/  H  i  n  \  i  i  l  l  ea/ 


2  r  i  /  \  c(k),cJi  i 

*  2^\bi  It  Pi  -  (ai  -  ac,i)  — "i 
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where 


w.  =-h.  + 

l  l 


V  0. 
n  i 


V  o.  tan  A 
n  l  ea 


+  b.  (  ■  J*l-m  +  a  .  -  a.  )  (k.  *  V  T.  tan  A  \ 
i  \  2n  c,i  1/  \  l  n  i  ea/ 


These  expressions  were  formulated  on  the  basis  of  hydrofoil  sections 
oriented  normal  to  the  hydrofoil  elastic  axis,  although  an  analogous 
procedure  would  apply  for  streamwise  sections.  If  the  modification  factor 
p^  were  set  equal  to  1  in  the  previous  expressions,  the  loading  would  be 
that  used  by  Yztes.  Three-dimensional  flow  effects  reduce  sectional  values 
of  p.  to  less  than  1. 


NONCIRCULATORY  MODIFICATION  FA'  -n 

In  attempting  to  devek  in  expression  for  the  spanwise  modification 
factor  p^,  the  literature  was  insulted  for  solutions  to  the  problem  of  a 
nonlifting  three-dimensional  body,  oscillating  in  a  fluid  medium  with  or 
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without  a  free  surface.  References  14  and  15  gave  some  general  approaches 

in  this  area;  however,  no  useful  analytical  expressions  have  been  derived. 

Therefore,  in  view  of  the  presence  of  virtual  mass  and  moment  in  the  non- 

circulatory  loading,  the  modification  factor  in  the  absence  of  a  free 

surface  was  inferred  from  the  distribution  of  virtual  mass  on  a  flat  plate 

oscillating  normal  to  its  own  plane.  The  factor  consisted  of  the  three- 

dimensional  virtual  mass  of  the  plate,  expressed  as  a  fraction  of  the  two- 

dimensional  virtual  mass,  outboard  of  the  spanwise  position  being  considered. 

This  quantity  may  be  expressed  as  a  spanwise  function  of  the  local  aspect- 

ratio  A.*  of  the  ith  section  in  the  form 
1 


A.* 


For  the  limiting  case  of  infinite  aspect  ratio,  corresponding  to  a  two- 

dimensional  hydrofoil,  p^  is  equal  to  1. 

A  different  modification  factor  is  required  for  high-frequency 

motion  when  a  free  surface  is  present.  By  referring  to  a  solution  of  the 

two-dimensional  problem  of  a  cylinder  swaying  in  a  fluid  medium  while 

piercing  the  free  surface,16  it  was  found  that  the  pressure  loading  on  a 

2 

thin  plate  was  strongly  dependent  on  the  nondimensional  frequency  (ta  /g. 

2 

As  shown  in  Figure  3,  large  values  of  &uj  / g  result  in  pressure  distribu¬ 
tions  which  decrease  to  zero  at  the  free  surface,  while  the  pressure 

approaches  the  free  surface  near  its  maximum  value  for  small  values  of 
2  2 

£oj  /g.  Therefore  for  values  of  lu  /g  of  10  or  greater,  a  free  surface 

effect  must  be  included  in  the  modification  factor.  The  free  surface  may 

2 

be  considered  a  reflecting  plane  for  values  of  fcw  /g  of  1  or  less.  For 

2 

hydrofoils  in  general,  Hoj  /g  is  very  large;  hence,  the  distribution  of  p. 

2  1 

should  be  similar  to  the  curve  with  Sxo  /g  =  00 .  To  simulate  this  distribu¬ 
tion  in  the  present  report,  Equation  (9)  was  used  to  calculate  the  values 
of  by  considering  the  midspan  of  the  hydrofoil  to  be  a  reflecting  plane. 

Total  virtual  mass  values  given  by  the  low-frequency  version  of  the 
formula  agreed  well  with  measured  values  reported  in  Reference  8,  while 
the  high-frequency  version  showed  only  fair  agreement.  This  comparison 
is  shown  in  Figure  4.  Nc  frequency  data  were  given  in  Reference  8  to 
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enable  determination  of  the  experimental  nondimens ional  frequency.  Average 
values  of  the  function  p^  over  each  hydrofoil  section  were  used  in  the 
flutter  computation. 

Figure  S  shows  two  distributions  of  p^.  The  curve  shown  for  the 
SwRI  30-inch  model  corresponds  to  a  reflecting  plane  at  the  root  of  the 
model.  The  curve  shown  for  Grumman  Model  1  has  incorporated  a  free-surface 

condition  by  considering  the  midspan  of  the  foil  to  be  a  reflecting  plane. 

2 

The  value  of  Au  /g  for  Grumman  Model  1  was  32  at  flutter,  and  therefore 

2 

the  free  surface  condition  was  appropriate.  When  fcoi  /g  is  not  known  for  a 
surface-piercing  model  at  flutter,  a  value  may  be  assumed,  using  perhaps 
the  first  torsional  frequency.  If  the  calculated  flutter  frequency  were 
significantly  different,  the  free-surface  boundary  condition  could  be 
changed  and  the  flutter  characteristics  recalculated. 

STEADY  LIFT  AND  AERODYNAMIC  CENTER  DISTRIBUTIONS 

In  accordance  with  the  procedure  used  by  Yates, ^  the  spanwise 
distributions  of  lift  and  aerodynamic  center  of  the  hydrofoils  were  calcu¬ 
lated  by  lifting-surface  theory.  Although  unsteady  lift  and  aerodynamic 
center  can  be  calculated  by  lifting-surface  theory,  only  steady  loading 
distributions  can  be  used  in  the  Yates  unsteady  loading  formulation.  There¬ 
fore  the  results  obtained  are  a  quasi-steady  approximation  to  the  exact 
unsteady  loading  characteristics. 

A  lifting-surface  theory  developed  and  programmed  for  the  IBM  7090 
computer  by  Widnall^  has  been  used  in  the  present  work.  This  theory  can 
treat  three-dimensional  planar  and  nonplanar  foil  configurations  that  are 
assumed  to  be  of  zero  thickness  and  at  small  angles  of  attack.  Surface¬ 
piercing  and  finite-depth  configurations  at  infinite  Froude  number  can  be 
treated  as  well  as  infinite-depth  configurations. 

The  Widnall  theory  has  not  been  fully  evaluated  for  agreement  with 
either  other  calculations  or  experiments  or  for  numerical  convergence. 
However,  available  calculrt'ons  agree  with  other  lifting-surface  theories 
and  are  accurate  at  low-asj.  ct  ratios  (less  than  3)  and  inaccurate  at  high- 
aspect  ratios  (more  than  3)  when  compared  to  experiments.  Some  numerical 
convergence  difficulties  also  occur  at  high  aspect  ratios.  In  view  of 
the  agreement  with  other  theories  and  the  extreme  versatility  of  the  Widnall 


theory,  it  has  been  used  for  the  present  steady  loading  calculations.  It 
is  recommended  that  the  sensitivity  of  the  flutter  calculation  to  variations 
in  lift  slope  and  center  of  pressure  be  determined,  particularly  at  high 
aspect  ratio. 

Steady  loading  was  required  for  each  flutter  model  in  the  form  of 
lift  slope  and  aerodynamic  center  as  functions  of  spanwise  position.  These 
quantities  were  calculated  by  assuming  a  uniform  angle  of  attack  along  the 
entire  span  of  the  foil.  While  it  would  be  possible  to  simulate  a  calculated 
or  measured-mode  shape  by  choosing  a  varying  angle  of  attack,  the  lift  slope 
and  aerodynamic  center  for  a  first  torsion  mode  were  similar  to  those  for 
the  uniform  angle  over  most  of  the  foil.  Effects  on  loading  due  to  varying 
deflection  were  therefore  regarded  as  second  order  effects  suitable  for 
future  investigation. 

Figure  6  shows  the  load  characteristics  calculated  for  the  five 

flutter  models-  All  models  except  the  SwRI  30-inch  model  were  only  partially 

submerged,  so  their  load  distributions  began  at  the  spanwise  position 

corresponding  to  the  waterline.  The  boundary  condition  at  the  free  surface 

was  determined  in  the  same  manner  as  for  the  noncirculatory  load-modification 

factor  described  previously.  For  values  of  the  nondimens ional  frequency 
2 

dm  /g  much  greater  than  1,  a  high-frequency  condition,  the  antisymmetrical 

boundary  condition  was  invoked,  resulting  in  lift  slopes  that  dropped  to 

2 

zero  at  the  free  surface.  For  2w  /g«l,  or  low  frequencies,  the  free  sur¬ 
face  was  treated  as  a  reflecting  plane  which  supported  symmetrical  loading. 
Ultimately,  the  frequency  used  for  this  determination  would  be  the  calculated 
flutter  frequency.  As  described  previously,  a  preliminary  estimate  could 
be  made  using  the  first  torsional  frequency,  subject  to  the  outcome  of  the 
flutter  calculation.  It  was  determined  that  Grumman  Models  1,  2,  and  3 
had  high-frequency  load  characteristics,  while  Grumman  Model  A  had  low- 
frequency  characteristics.  The  SwRI  30-inch  model  was  attached  to  a  root- 
reflecting  plate  during  testing. 

STRUCTURAL  VIBRATION  CHARACTERISTICS 

Natural  frequencies  and  associated  mode  shapes  of  the  five  models 
selected  in  this  study  were  computed.  The  mode  shapes  were  expressed  in 
terms  of  nodal  lines.  Both  analytical  results  and  experimental  data  are 
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presented  in  the  following  sections.  The  analytical  calculations  show  very 
good  agreement  with  experimental  results.  Hence,  it  is  felt  that  the 
mathematical  model  used  for  the  flutter  analysis  adquately  represents  the 
structural  characteristics  of  the  models. 

SwRI  30- INCH  FLUTTER  MODEL 

Table  1  gives  the  calculated  and  measured  natural  frequencies  for 
the  SwRI  30-inch  flutter  model  in  vacuum,  and  Table  2  gives  the  calculated 
and  measured  natural  frequencies  in  water.  Table  3  shows  the  calculated 
and  measured  first  coupled  bending  mode  shape  in  vacuum  and  first  coupled 
torsional  mode  shape  in  vacuum.  Figures  7  and  8  show  the  calculated  nodal 
lines  in  vacuum  and  in  water,  respectively.  No  experimental  results  on 
measured  nodal  lines  are  available. 

GRUMMAN  MODEL  1 

Table  4  gives  the  calculated  and  measured  natural  frequencies  for 
Grumman  Model  1  in  vacuum.  Table  5  gives  the  analytical  natural  frequencies 
in  water;  no  test  data  are  available.  For  this  particular  model,  the 
elastic  axis  coincides  with  the  center  of  gravity,  and  is  also  at  the 
midchord  of  the  hydrofoil.  Therefore,  at  zero  speed,  it  is  a  decoupled 
bending  and  torsion  system,  either  in  vacuum  or  in  water.  Figures  9  and  10 
show  the  calculated  nodal  lines  in  vacuum  and  in  water,  respectively.  Only 

g 

the  nodal  lines  in  air  have  been  measured,  and  they  are  included  in 
Figure  9. 

GRUMMAN  MODEL  2 

Table  6  gives  calculated  and  measured  natural  frequencies  for 

Grumman  Model  2  in  vacuum.  Table  7  gives  the  analytical  natural  frequencies 

in  water;  no  test  data  are  available.  Figures  11  anc'  12  give  the  analytical 

nodal  lines  in  vacuum  and  in  water,  respectively.  Only  the  nodal  lines  in 

8 

air  have  been  measured,  and  they  are  included  in  Figure  11.  It  is  inter¬ 
esting  to  note  that  the  order  of  the  fifth  and  sixth  modes  in  vacuum  was 
reversed  when  the  model  was  in  water. 
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GRUMMAN  MODEL  3 


Tables  8  and  9  give  the  analytical  natural  frequencies  for  Grumman 
Model  3  in  vacuum  and  in  water,  respectively.  The  analytical  nodal  lines 
in  vacuum  and  in  water  are  also  shown  in  Figures  13  and  14,  respectively. 
No  experimental  data  are  available  for  either  natural  frequencies  or  nodal 
1 ines . 

GRUMMAN  MODEL  A 

Calculated  and  measured  natural  frequencies  for  Grumman  Model  A 
in  vacuum  are  given  in  Table  10.  Analytical  natural  frequencies  in  water 
are  given  in  Table  11;  no  test  data  are  available.  The  calculated  nodal 
lines  in  vacuum  and  in  water,  along  with  the  measured  nodal  lines  in  air, 
are  shown  in  Figures  15  and  16,  respectively. 

FLUTTER  RESULTS 


SwRI  30- INCH  FLUTTER  MODEL 

The  SwRI  30- inch  flutter  model  underwent  structural  failure  due  to 
flutter  at  a  speed  of  48.1  knots  and  a  frequency  of  17.5  Hz,  as  described 
in  Reference  4.  Measured  values  of  exponential  decay  factor  and  oscillation 
frequency  for  the  unstable  mode  are  shown  in  Figures  17  and  18,  respectively. 

Many  flutter  calculations  have  been  published  for  this  model.  Fig¬ 
ures  17  and  18  show  the  results  of  four  published  calculations  as  Theories  1 
through  4.  These  calculations  contained  no  empirical  corrections.  All 
four  of  the  theories  gave  drastically  unconservative  predictions  of  flutter 
speed,  the  closest  agreement  with  experiment  being  obtained  from  the  gener¬ 
alized  lifting  surface  theory  of  Reference  3  at  67.5  knots. 

In  contrast  to  the  unconservative  predictions  of  the  previously 
described  theories,  two  calculations  based  on  Yates  hydrodynamic  loading 
gave  conservative  results.  The  calculation  using  unmodified  Yates  loading, 
identified  as  Theory  5  in  Figures  17  and  18,  predicted  an  overconservative 
flutter  speed  of  30  knots.  It  is  felt  that  the  difference  between  the 
present  result,  which  was  conservative  by  38  percent,  and  the  original 
Yates  result, ^  which  was  conservative  by  20  percent,  was  due  to  numerical 
approximations  accentuated  by  the  nearly  horizontal  orientation  of  the 
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damping  curve.  An  improved  result,  shown  as  Theory  6,  was  obtained  using 
modified  Yates  loading.  The  modified  Yates  theory  gave  a  flutter  speed 
of  40. S  knots,  which  was  16  percent  conservative,  and,  in  contrast  to  all 
of  the  other  theories,  showed  good  agreement  with  measured  damping  values. 
Both  of  these  calculations  assumed  a  value  of  0.02  for  the  structural 
damping  parameter  g.  All  of  the  predicted  frequencies  followed  the  measured 
frequencies  fairly  well  over  the  observed  range,  and  corresponded  to  the 
second  mode  of  response. 

GRUMMAN  MODELS  1,  2,  AND  3 

A  description  of  the  flutter  testing  of  Grumman  Models  1,  2,  and  3 
is  given  in  Reference  8.  The  models  were  accelerated  until  flatter 
occurred;  the  speed  at  which  flutter  stopped  during  deceleration  was  taken 
as  the  flutter  speed.  Figures  19  through  21  show  the  observed  flutter 
characteristics. 

Three  published  flutter  speed  predictions  for  Model  1,  shown  as 
Theories  1,  2,  and  3  in  Figure  19,  were  grossly  unconservative  when  com¬ 
pared  with  the  experimental  flutter  speed  of  66  knots.  Both  of  the  flutter 
calculations  based  on  Yates  hydrodynamic  loading,  shown  as  Theories  4  and  5 
in  Figure  19,  were  similarly  unconservative.  Using  unmodified  Yates  load¬ 
ing,  a  flutter  speed  of  118  knots  was  obtained,  which  was  79  percent  above 
the  experimental  value.  A  more  unconservative  result  was  obtained  with 
modified  Yates  loading  which  gave  a  flutter  speed  of  148  knots,  124  per¬ 
cent  above  experiment.  The  predicted  and  observed  frequencies  agreed 
fairly  well  as  functions  of  speed,  including  a  "dead  beat"  point  at  36 
knots,  obtained  from  the  test  log.*  Flutter  was  predicted  to  ~ccur  in  the 
first  mode  of  response. 

The  model  response  characteristics  calculated  for  Model  2  using  both 
Yates  and  modified  Yates  loading,  are  shown  in  Figure  20.  Flutter  occurred 
experimentally  at  81  knots.  The  Yates  analysis  predicted  a  flutter  speed 
of  87  knots,  which  was  unconservative  by  7  percent.  After  modification 


♦Private  communication  with  Mr.  Chailes  E.  Squires,  Jr.,  Grumman 
Aerospace  Corporation. 
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to  the  loading,  a  flutter  speed  of  104  knots  was  calculated,  28  percent 
above  experiment.  Frequency  predictions  for  Model  2  followed  the  observed 
frequencies  fairly  well  as  a  function  of  speed.  The  frequencies  between 
84  and  100  knots  were  obtained  from  the  experimental  log  book.*  Flutter 
was  again  piedicted  to  occur  in  the  first  mode  of  response. 

There  was  little  difference  in  the  response  characteristics  calcu¬ 
lated  for  Model  3  using  the  two  versions  of  Yates  loading  (see  Figure  21). 
The  calculations  gave  a  flutter  speed  of  61  knots  for  the  Yates  analysis 
and  a  slightly  higher  flutter  speed  of  63  knots  for  the  modified  Yates 
analysis.  These  predictions  were  approximately  30  percent  above  the 
experimental  value  of  48  knots.  The  predicted  frequency  curves  agree 
fairly  well  with  the  frequency  observed  during  flutter. 

GRUMMAN  MODEL  A 

Figure  22  shows  flutter  test  results  for  Grumman  Model  A  along  with 
two  previously  published  flutter  predictions  and  the  two  present  Yates 
analyses.  The  flutter  results  were  taken  from  Reference  9  and  corresponded 
to  a  sweep  angle  of  5  degrees. 

Theory  1  in  Figure  22,  based  on  two-dimensional  Theodorsen  loading, 
agreed  exactly  with  the  experimental  flutter  speed  of  32.5  knots.  Theory  2, 
a  generalized  lifting  surface  theory,  gave  an  unconservative  flutter  speed 
of  47.2  knots,  45  percent  above  experiment.  The  Yates  analysis,  shown  as 
Theory  3,  was  14  percent  unconservative  in  its  prediction  of  37  knots.  A 
low-frequency  free  surface  condition  was  assumed  for  this  calculation. 

The  modified  Yates  analysis  was  made  with  two  free  surface  condi¬ 
tions  for  comparison.  As  previously  described,  a  low  frequency  surface¬ 
piercing  foil  should  be  treated  as  though  the  free  surface  were  a  reflecting 

2 

plane.  Because  the  nondimensional  frequency  Hw  /g  was  0.35  for  Model  A  at 
flutter,  a  low-frequency  condition  was  assumed  in  calculating  the  theoret¬ 
ical  Curve  4a  in  Figure  22.  An  unconservative  flutter  speed  of  38  knots 
was  obtained,  17  percent  above  the  experimental  value.  When  the  high- 
frequency  surface  condition  was  used  in  calculating  Theory  4b,  the  less 


*Refer  to  footnote  on  page  13. 
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accurate  flutter  speed  of  42  knots  obtained  was  29  percent  unconservative. 

It  is  concluded  that  the  free  surface  condition  to  be  used  in  calculating 
the  noncirculatory  loading  modification  should  be  determined  by  the  non- 
dimensional  frequency  of  the  hydrofoil  at  flutter. 

All  theoretical  approaches  except  Theory  2  gave  accurate  frequency 
predictions  at  the  speed  at  which  flutter  occurred  experimentally.  As 
obtained  for  the  other  surface-piercing  struts,  the  analyses  based  on 
different  forms  of  Yates  loading  predicted  that  flutter  instability  would 
occur  in  the  first  mode  of  response. 

DISCUSSION 

The  two  Yates  flutter  analyses  did  not  show  general  agreement  with 
experimental  results.  Conservative  predictions  were  obtained  for  the 
SwRI  30-inch  model,  a  fully-wet;  d  strut  with  a  mass  ratio  of  0.99.  Uncon¬ 
servative  predictions  were  obtained  for  Grumman  Models  1,  2,  3,  and  A; 
these  were  surface-piercing  struts  with  mass  ratios  between  0.09  and  0.28. 
The  modified  Yates  analysis  always  gave  a  higher  flutter  speed  than  the 
unmodified  analysis.  Therefore  the  noncirculatory  load  modification 
improved  the  results  for  the  SwRI  model  and  worsened  the  predictions  for 
the  surface-piercing  struts. 

An  interesting  pattern  was  observed  in  the  modes  of  response  of 
the  flutter  models.  The  conservative  and  unconservative  flutter  speed 
predictions  differed  in  the  predicted  flutter  mode.  When  conservative 
results  were  obtained,  flutter  was  predicted  to  occur  in  the  second  mode 
of  response,  with  deflections  primarily  resembling  the  fundamental  torsion 
mode  of  a  cantilever  beam.  Unconservative  results  were  always  associated 
with  the  first  mode  of  response,  resembling  the  fundamental  bending  mode 
of  a  cantilever  beam.  The  predicted  frequencies  agreed  well  with  observed 
frequencies,  and  therefore  it  is  felt  that  the  predicted  flutter  modes 
were  the  ones  that  actually  occurred.  The  tendency  of  the  two  Yates 
analyses  to  be  conservative  or  unconservative  in  flutter  predictions  appears 
to  be  associated  with  the  mode  of  response  at  flutter. 

Several  hydrodynamic  factors  differed  among  the  models,  but  they 
were  not  fel*  to  be  responsible  for  the  differences  in  the  flutter  modes. 
These  hydrodynamic  factors  were  mass  ratio,  submergence,  sweep  angle,  and 
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profile.  Rather  than  these  factors,  it  is  felt  that  differences  in 
structural  properties  of  the  models  caused  the  variation  in  flutter  modes. 
The  structural  differences  may  be  observed  in  the  two  groupings  of  struc¬ 
tural  modes  which  appear  in  Figures  7  through  16.  The  SwRI  30-inch  model 
displayed  a  first  bending  mode  followed  by  a  first  torsion  mode,  while  the 
Grumman  models  displayed  a  first  bending  mode  followed  by  a  second  bending 
mode.  The  latter  models  displayed  a  first  torsion  mode  in  the  third 
position.  These  two  mode  sequences  occurred  both  in  vacuum  and  in  water. 
The  mode  sequence  variation  was  caused  by  a  difference  in  the  relative 
stiffness  of  bending  and  torsion,  and  in  the  mass  distribution,  of  the  two 
types  of  models.  The  SwRI  30-inch  model  had  a  relatively  greater  tendency 
to  oscillate  in  torsion  than  the  other  models. 

On  the  basis  of  the  above  correlation,  it  is  concluded  that  the 
Yates  analysis  and  the  modified  Yates  analysis  may  be  useful  in  obtaining 
conservative  hydrofoil  flutter  speed  predictions  for  flutter  which  occurs 
in  a  torsional  mode.  Further  correlation  with  data  is  required  to  confirm 
this  conclusion.  It  does  not  appear  that  these  analyses  are  useful  for 
predominantly  bending  flutter.  Additional  investigation  of  the  hydro¬ 
dynamics  is  required  in  order  to  obtain  improved  results  for  bending 
flutter.  In  particular,  this  investigation  should  include  a  determination 
of  the  sensitivity  of  the  flutter  calculation  to  variations  in  spanwise 
input  parameters.  Important  parameters  should  then  be  checked  with  avail¬ 
able  data.  The  entire  hydrodynamic  loading  formulation  should  be  compared 
with  data  given  in  References  18  and  19. 

None  of  the  other  theories  presented  except  the  two-dimensional 
Theodorsen  gave  conservative  predictions  for  any  of  the  hydrofoils.  The 
Theodorsen  theory  gave  exact  agreement  for  a  bending  flutter  occurrence  but 
an  unconservative  prediction  for  the  torsional  case.  Therefore  none  of  the 
other  theories  appear  to  be  useful  in  treating  hydrofoil  flutter  of  either 
the  bending  or  the  torsional  type. 

CONCLUSIONS 

1.  The  Yates  analysis  and  the  modified  Yates  analysis  may  be  useful 
in  obtaining  conservative  hydrofoil  flutter  speed  predictions  for  flutter 
which  occurs  in  a  torsional  mode.  It  does  not  appear  that  these  analyses 

are  useful  for  predominantly  bending  flutter. 
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2.  Both  Yates  analyses  are  capable  of  giving  fairly  accurate 
predictions  of  frequency  as  a  function  of  speed. 

3.  The  other  theories  used  for  comparison  with  the  Yates  theories 
do  not  appear  to  be  useful  in  treating  hydrofoil  flutter  of  either  the 
bending  or  the  torsional  types. 

4.  Structural  characteristics  of  the  flutter  models  were  adequately 
represented  by  the  lumped -parameter  mathematical  model. 

5.  The  free-surface  condition  to  be  applied  in  calculating  the 
noncirculatory  loading  modification  for  a  surface-piercing  hydrofoil  should 
be  determined  according  to  the  nondimensional  frequency  of  the  hydrofoil 

at  flutter. 


RECOMMENDATIONS 

1.  Further  correlation  with  data  is  recommended  to  establish  the 
validity  of  the  Yates  theories  in  treating  predominantly  torsional  flutter, 
as  compared  to  predominantly  bending  flutter,  of  hydrofoils. 

2.  Additional  investigation  of  the  hydrodynamics  is  required  in 
order  to  obtain  improved  resuits  for  bending  flutter.  In  particular,  this 
investigation  should  include  a  determination  of  the  sensitivity  of  the 
flutter  calculation  to  variation  in  spanwise  input  parameters.  Important 
parameters  should  then  be  checked  with  available  data.  The  entire  hydro- 
dynamic  loading  formulation  should  be  compared  with  the  data  given  in 
References  18  and  19. 
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Figure  6b  -  Spanwise  Distribution  cf  Lift  Coefficient  Slope 
Figure  6  -  Steady  Hydrodynamic  Loading  Characteristics 
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Figure  13  -  Analytical  Nodal  Lines  in  Vacuum  for  Grumman  Model  3 
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Figure  14  -  Analytical  Nodal  Lines  in  Water  for  Grumman  Model  3 
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Figure  15  -  Analytical 
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Figure  19b  -  Frequency- 

Figure  19  -  Critical  Flutter  Modes  as  Functions 
of  Speed  for  Grumman  Model  1 
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FREQUENCY  IN  HERTZ 


Figure  20  -  Critical  Flutter  Modes  as  Functions 
of  Speed  for  Grumman  Model  2 
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Figure  21b  -  Frequency 


Figure  21  -  Critical  Flutter  Modes  as  Functions 
of  Speed  for  Grumman  Model  3 
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Figure  22b  -  Frequency 

Figure  22  -  Critical  Flutter  Modes  as  Functions 
of  Speed  for  Gruiranan  Model  A 
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TABLE  1 

Comparison  of  Analytical  and  Measured  Natural  Frequencies 
in  Vacuum  of  SwRI  30-Inch  Flutter  Model 


Mode  Number 

4 

Measured  in  Air 

Hz 

Analytical  in  Vacuum 

Hz 

1 

9.63 

10.1 

2 

31.9 

28.3 

3 

50.2 

4 

79.4 

5 

116.7 

6 

126.4 

TABLE  2 

Comparison  of  Analytical  and  Measured  Natural 
Frequencies  in  Water  of  SwRI  30- Inch 
Flutter  Model 

4  | 

u  .  ..  .  Measured  Analytical 

Mode  Number  Hz  Hz 


TABLE  3 


Comparison  of  Analytical  and  Measured  First  Bending  and  First 
Torsional  Mode  Shapes  in  Vacuum  of  SwRI  30-Inch  Flutter  Model 


y/L 

First  Coupled  Bending-Node  Shape 
Translational  Deflection 

First  Coupled  Torsional-Node  Shape 
Rotational  Deflection 

Measured 
in  Air* 

Analytical 

in  Vacuum 

Measured 
in  Air* 

Analytical 
in  Vacuum 

Root 

0 

mm 

0. 

0. 

0. 

0.1 

0.018 

0.213 

0.202 

0.2 

msm 

0.065 

0.403 

0.393 

0.3 

0.139 

0.570 

0.564 

0.4 

1  ESI  1 

0.232 

0.710 

0.708 

0.5 

0.342 

0.825 

0.823 

0.6 

0.488 

0.464 

0.915 

0.905 

0.610 

0.593 

0.966 

0.959 

0.8 

0.738 

0.727 

0.990 

0.988 

0.870 

0.863 

0.999 

0.999 

Tip 

1.0 

1.000 

1.000 

1.000 

TABLE  4 

Comparison  of  Analytical  and  Measured  Natural 
Frequencies  in  Vacuum  of  Grumman  Model  1 


Measured  in  Air 
Hz 

Analytical  in  Vacuum 
Hz 

Mode  Shape 

Ref.  8 

Ref.  3 

i 

1 

2.87 

2.83 

2.88 

1st  Bending 

2 

18.2 

18.00 

17.43 

2nd  Bending 

3 

37.5 

38.08 

38.59 

1st  Torsion 

4 

49.7 

49.92 

46.70 

3rd  Bending 

5 

94.8 

97.02 

86.23 

4th  Bending 

6 

111.0 

116.0 

114.8 

2nd  Torsion 
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TABLE  5 


Analytical  Natural  Frequencies  in  Water 
of  Grumman  Model  1 


Node  Number 

Analytical  In  Water 

H2 

Node  Shape 

1 

1.38 

1st  Bending 

2 

8.39 

2nd  Bending 

3 

21.14 

1st  Torsion 

4 

23.39 

3rd  Bending 

5 

44.59 

4th  Bending 

6 

66.00 

2nd  Torsion 

TABLE  6 

Comparison  of  Analytical  and  Measured  Natural  Frequencies 
in  Vacuum  of  Grumman  Model  2 


Node 

Number 

Q 

Neasured  in  Air 

Hz 

Analytical  in  Vacuum 
l.U 

Predominant 
Node  Shape 

1 

4.90 

4.82 

1st  Bending 

2 

79.8 

29.03 

2nd  Bending 

3 

47.5 

42.69 

1st  Torsion 

4 

80.82 

77.18 

3rd  Bending 

5 

138.0 

127.4 

2nd  Torsion 

6 

" 

140.6 

4th  Bending 

TABLE  7 

Analytical  Natural  Frequencies  in  Water 
of  Grumman  Model  2 


Node 

Numoer 

Analytical  in  Water 

Hz 

Predominant 
Node  Shape 

1 

2.58 

1st  Bending 

2 

15.61 

2nd  Bending 

3 

29.22 

1st  Torsion 

4 

43,72 

3rd  Bendinq 

5 

81.45 

4th  Bending 

6 

92.66 

2nd  Torsion 

TABLE  9 


Analytical  Natural  Frequencies  in  Water 
of  Grumman  Model  3 


TABLE  10 

Comparison  of  Analytical  and  Measured  Natural  Frequencies 
in  Vacuui.i  of  Grumman  Model  A 


Mode 

Number 


Measured  in  flir^ 

Hz 

Analytical  in  Vacuum 

Hz 

Predominant 
Mode  Shape 

3.20 

3.21 

1st  Bending 

19.67 

19.31 

2nd  Bending 

32.45 

33.90 

1st  Torsion 

53.88 

51.22 

3rd  Bending 

107.0 

92.85 

4th  Bending 

100.7 

101.5 

2nd  Torsion 

-- 

140.2 

5th  Bending 

164.9 

166.9 

3rd  Torsion 
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TABLE  11 

Analytical  Natural  Frequencies  in  Water 
of  Grumman  Model  A 


MGde 

Number 

Analytical  in  Water 

Hz 

Predominent 
Mode  Shape 

1 

1.11 

1st  Bending 

2 

6.47 

2nd  Bending 

3 

14.46 

1st  Torsion 

4 

17.19 

3rd  Bending 

5 

31.80 

4th  Bending 

6 

42.95 

2nd  Torsion 

7 

49.24 

5th  Bending 

8 

67.60 

_ 

6th  Bending 
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APPENDIX 


DESCRIPTION  OF  FLUTTER  MODELS 


SwRI  30- INCH  MODEL 

A  detailed  description  of  this  model  was  given  in  Reference  4.  The 
model  consisted  of  a  steel  spar  to  which  lead  segments  were  attached.  The 
planform  and  profile  are  shown  in  Figure  23,  and  the  structural  character¬ 
istics  are  summarized  in  Table  12.  A  reflecting  plate  was  attached  to  the 
root  of  the  model  to  eliminate  the  free  surface  effect. 

GRUKNAN  MODEL  1 

This  model  was  Grumman  Model  1  of  Reference  8.  It  was  constructed 
of  solid  steel.  The  planform  and  profile  are  shown  in  Figure  24,  and  the 
structural  characteristics  are  summarized  in  Table  13. 

GRUMMAN  MODELS  2  AND  3 

These  two  models  were  Grumman  Models  2  and  3  of  Reference  8. 

Models  2  and  3  (Figure  25)  are  geometrically  identical  struts  constructed 
of  solid  steel  and  solid  aluminum  respectively.  Tneir  structural  charac¬ 
teristics  are  summarized  in  Tables  14  and  15. 

GRUMMAN  MODEL  A 

This  model  was  Grumman  Model  A  of  Reference  9.  The  planform  and 
profile  are  shown  in  Figure  26  and  its  structural  characteristics  are 
summarized  in  Table  16. 


Preceding  page  blank 
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Grumman  Models  2  and  3  Figure  26  -  Grumman  Model  A 


TABLE  12 

SwRI  30-Inch  Flutter  Model  Structural  Characteristics 


Model  Parameter 

Measured  Value 

Aspect  ratio  AR 

5.0 

Semi  span  L  in  in. 

30.0 

Sweep  angle  Aga  in  deg 

0.0 

Semi  chord  b  in  in. 

6.0 

2  2 

Mass  per  unit  span  m  in  lb-sec  /in. 

0.0105 

Mass  ratio  u 

0.99 

Elastic  axis  location  a 

-0.5 

Center  of  gravity  location  xa 

0.524 

2 

Radius  of  gyration  r 

a 

0.512 

2 

Bending  stiffness  El  in  lb-in. 

3.40  x  106 

2 

Torsional  stiffness  GJ  in  lb-in. 

0.973  x  106 

Total  weight  (wing  only)  in  lb 

121.2 
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TABLE  13 


Grumman  Model  1  Structural  Characteristics 


Model  Parameter 


Value 


Aspect  ratio  AR 
Semispan  L  in  in. 
Sweep  anqle  A  in 

6  n 

Semichord  b  in  in. 


10.35 

62.1 


deq 


15.0 


6.0 


2  2 

Mass  per  unit  span  m  in  lb-sec  /in. 

Mass  ratio  u 

Elastic  axis  location  a 


0.00232 

0.22 

0.0 


Center  of  qravity  location  x 

01 


o.o 


Radius  of  qyration  r 


0.1975 


p 

Bendinq  stiffness  El*  in  lb-in. 
Torsional  stiffness  GJ*  in  lb-in. 


0.9564  x  106 
1.555  x  106 


*E!  is  0.928  x  106  lh-in. Z.  GJ  is  1.52  x  106  lb-in.2  in 
Reference  3,  in  which  results  were  obtained  by  calibratinq 
after-flutter-tested  model. 


TABLE  14 

Grumman  Model  2  Structural  Characteristics 


Model  Parameter 

Value 

Aspect  ratio  AR 

8.57 

Semi  span  l  in  in. 

49.693 

Sweep  anqle  A„,  in  deq 

6o 

15 

Semi  chord  b  in  in. 

5.796 

2  2 

Mass  per  unit  span  m  in  lb-sec  /in. 

0.00273 

Mass  ratio  u 

0.28 

Elastic  axis  location  a 

0.38 

Center  of  qravity  location  x^ 

-0.214 

2 

Radius  of  qyration  rft 

0.358 

9 

Bendinq  stiffness  El  in  lb-in." 

1.293  x  106 

2 

Torsional  stiffress  Gd  in  lb-in. 

2.070  x  106 
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TABLE  15 


Grumman  Model  3  Structural  Characteristics 


Model  Parameter 

Value 

Aspect  ratio  AR 

8.57 

Semispan  L  In  In. 

49.693 

Sweep  anqle  A,,,  In  deg 

6a 

15 

Semichord  b  in  in. 

5.796 

2  2 

Mass  per  unit  span  m  in  lb-sec  /in. 

0.00094 

Mass  ratio  y 

0.10 

Elastic  axis  location  a 

0.38 

Center  of  gravity  location  x^ 

-0.214 

2 

Radius  of  gyration  r 

<x 

0.326 

2 

Bending  stiffness  El  in  lb-in. 

0.431  x  10f 

2  i 

Torsional  stiffness  GJ  in  lb-in. 

0.691  x  106 

TABLE  16 

Grumman  Model  A  Structural  Characteristics 

Model  Parameter 

Value 

Aspect  ratio  AR 

8.251 

Semi  span  L  in  in. 

48 

Sweep  angle  A  in  deg 

5 

Semi  chord  b  in  in. 

6 

2  .  2 

Mass  per  unit  span  in  in  lb-sec  /in. 

0.00097475 

Mass  ratio  y 

0.09 

Elastic  axis  location  a 

-0.0784 

Center  of  gravity  location  xft 

0.2724 

2 

Radius  of  gyration  rft 

0.3254 

2 

Bending  stiffness  El  in  lb- in. 

0.174  x  106 

2 

Torsional  stiffness  GJ  in  lb-in. 

0.370  x  106 
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